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Abstract 

We prove the equivalence of ensembles for Bernoulli measures on Z conditioned 
on two conserved quantities under the situation that one of them is spatially inho- 
mogeneous. For the proof, we extend the classical local limit theorem for a sum of 
Bernoulli independent sequences to those multiplied by linearly growing weights. The 
motivation comes from the study of random Young diagrams. We discuss the relation 
between our result and the so-called Vershik curve which appears in a scaling limit for 



height functions of two-dimensional Young diagrams. 



1 Introduction 
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The equivalence of ensembles, that is the asymptotic equivalence of canonical and grand 
canonical ensembles for large systems, plays a fundamental role in equilibrium statistical 
physics [13], [6], [9] and also in some problems related to statistics [2], [3]. It is mostly 
discussed for canonical ensembles obtained under conditioning on spatially homogeneous 
physical quantities. In our problem, the grand canonical ensembles are simply Bernoulli 
measures on Z, but the system has two conservation laws so that the corresponding canon- 
ical ensembles are defined through conditioning on a quantity which is not translation- 
invariant. As a result, the macroscopic profile for the grand canonical ensemble turns out 
to be spatially dependent. We will show that this profile has a connection to the so-called 
Vershik curve which appears in a scaling limit for two-dimensional Young diagrams, |14| . 

ID. 01- 

The height function ip q of a two-dimensional Young diagram, which is associated with 
a distinct partition q = {qi > q2 > • • • > qx > 1} of a positive integer M by positive 
integers {qi\f = i (i.e. M = ^i=i%); i s given by a right continuous non-increasing step 
function 

K 



Graduate School of Mathematical Sciences, The University of Tokyo, Komaba, Tokyo 153-8914, Japan. 

e-mail: funaki@ms.u-tokyo.ac.jp, Fax: +81-3-5465-7011. 
Keywords: equivalence of ensembles, local limit theorem, Young diagram, Vershik curve. 
Abbreviated title [running head): Equivalence of ensembles and applications to random Young diagrams 
2010MSC: primary 60K35, secondary 82B20, 60D05, 60F05. 
Supported in part by the JSPS Grants (A) 22244007 and 21654021. 



1 



Its height difference rj = {ryfc}fc £ N,N = {1, 2, . . .}, is defined by 



Vk = - 1) - ip q (k) e {0,1} or rj k = 1 



fe 6 {q l } 



K 

i=l- 



Then, in terms of rj, the height K of the Young diagram at u = is represented as 



The sequence rj = {rj^keN can be interpreted as defining a configuration of particles on 
N in the sense that the site k is occupied by a particle if rjk = 1 and vacant if rjk = 
and, from the viewpoint of random Young diagrams, it is natural to restrict the space of 
configurations rj to those satisfying X^fceN Vk = K and X^fceN ^Vk = M, given two constants 
K and M. The first condition is equivalent to assign the total number of particles in the 
system, while the second sum involves a non-translation-invariant linearly growing weight 
k. The simplest random structure can be introduced to the set of Young diagrams with 
height K at u = and area M by means of a uniform probability measure on the space of 
configurations ij with these two constraints. In fact, such random structure is the subject 
of our study. 

The original motivation of this paper comes from the study of the hydrodynamic 
limit for an area-preserving random dynamic on two-dimensional Young diagrams with 
height differences restricted to be or 1. Such dynamic, which is a kind of a surface 
diffusion model studied in physics [12], is introduced in [5]. As we have seen above, it can 
be transformed into an equivalent particle system on N with two conservation laws. 

We formulate our problem and main results in Section 2. For the proof, we need to 
extend the local limit theorem for a sum of independent random variables. If the random 
variables satisfy certain proper moment conditions (see, e.g., Theorems VII. 4 and VII. 12 
of [H]), the classical theory applies, but in our case, such conditions do not hold since the 
random variables have unbounded linearly growing weights. This is discussed in Section 3 
relying on the fact that the weights are linear by analyzing the behavior of characteristic 
functions. Sections 4 and 5 are devoted to the proofs of main results. We finally show 
that the limit curves obtained in [T] and here are identical, although the random structures 
imposed on the set of Young diagrams are different. 

2 Main results 

Let us state our problem precisely. We consider a particle system on Z rather than on 
N following the usual setting in statistical physics. Each site is occupied by at most one 
particle. Therefore the particle configuration on Z is represented by rj = {rjk}kez G S := 
{0, 1} Z , where rjk = 1 means that the site k is occupied by a particle and rjk = means 
that k is vacant. For a finite set A in Z and rj £ S, we define 



Vv(o) = J2 r,k 



and the area M as 
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For given K, M and let v\ t! K,M be the uniform probability measure on the configu- 

ration space T,A e ,K,M = {v £ ^A £ ; K\ t {rj) = K, M\ e (rj) = M}, where = {k G Z; | fe| < ■£} 
and Sa = {0, 1} A . For a G (0, 1) and a finite set A in Z, let f A be the Bernoulli measure 
on S A with mean a, that is, ^ A (r?) = (1 - a )tt k&A ™=°} a ${ k ^Vk=i} f or eac h ^ g £ A . Then, 



the measure u^ e ,K,M is identical to the conditional measure of v£ e or i/£ with A; C A on 
^AaK.m: For £ G Y ..\,.i<.\h 



M) 



= "a(VA e = £\Z Ae ,K,M X £ A \ A J 

= ^ofa = (f,»7A\Aj|£A*Jir,Af x {*7a\a<}) 

independently of the choice of the outer conditions 7/A\A^i see also Lemma B~Tj below. We 
have denoted by rj\ t and r/A\A £ the configurations n restricted on these sets. 

We first observe the possible values of K = K\ e {rj) and M = M\ e (n) for rj G S^. It 
is easy to see that < K < 21 + 1 , while the values of M range as 

|M| < 1(K(2£ + 1) -K 2 ) 

under the condition that K\ £ (r]) = K. Indeed the extreme values of M given K are 
attained when all K particles are closely packed at the right or left most edges on A^. 
Accordingly, we have that 



K 



M 



{21 + 1) 




K 



21 + 1 



We fix p G N and assume that, for £ G N and 1 < j < p, sequences K = Kg, M = 
and kj = kj ; t are given and satisfy 

(2-1) lim-5— = ^(0,1), 

£-»oo LI + 1 



M 

/Too (2/Ti) 



(2.2) Urn — ; = m G (-v/2, w/2), 



£->oo 



(2.3) lim ^ = x, G (-1,1), 



respectively, with distinct limits {xj}^ =1 , where v = p(l — p). A function / on E is called 
local if it depends only on finitely many coordinates in n. The shift operators t% are defined 
by Tif(r]) = firiq) and (Tin)k = r\k+i for i,k G Z and n G S. We are now at the position 
to formulate our main theorem. 

Theorem 2.1. Let fj = fj(j]),l < j < p be local functions on S. Then, under the 
conditions (|2.ip ~ (|2.3p . we /iaue that 

p v 

3=1 3=1 
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where v a ,a G (0,1) denotes the Bernoulli measure on S with mean a and 

(2.4) P(x) = P(x;a,b)= bx f*° -, *G[-1,1], 

e 0:c a + (1 — a) 

with two parameters a G (0, 1) and 6 G E determined from p and m by the relations: 

(2.5) y /?(x)cte = 2p, J x/3(x)dx = 4m. 

T/ie convergence is uniform in (K, M) in the region that e < Kj (2t + 1) < 1 — e and 
M/(2£ + l) 2 G (-v/2 + e,v/2- e) for every e > 0. 

For every p and m, one can find a and b uniquely as is shown in Lemma 14.21 below. 
This theorem asserts that, as i — > oo under the canonical ensemble VK t ,K,M^ the limit 
distributions are asymptotically independent for microscopic regions which are macro- 
scopically separated, and the microscopic limit distribution around the macroscopic point 
x G (—1,1) is the grand canonical ensemble vpt x ) with macroscopically dependent profile 
/3(x). It will be useful to have another expression of /?: 

(2.6) j 9(x) = i^M = i(Iog»(x))' I 

where g(x) = e bx a + (1 — a). In particular, we have that 

1 e b a + (1 - a) 



f 1 1 p b a A 

(2.7) / ^(x)dx=ilog4F 

7-1 b e- b a 



(1-a)' 

1 #)i 4 log ( ebQ+ ( 1 - Q ^ + ( 1 - a » 

• (1 — a) z 



-i 



1 f. / ae b \ , / ae 



• 1-^ 

where ^(-z) := — Jq 7 1°§(1 — *)^*> ^ < 1 is the Euler dilogarithm, see [7j, p. 642. 

Remark 2.1. (1) Theorem \2.1\ gives the equivalence of ensembles under the situation that 
the system has two conservation laws, especially, one of them is not translation-invariant. 

(2) The uniformity of the convergence near the boundary values of K/(2£+l) and M/ (2£+ 
l) 2 can not be shown, because the cumulant \$, which controls the error estimate in the 
local limit theorem, diverges near the boundary values. 

(3) The macroscopic mean (3(x) of the limit measure is not translation-invariant, but the 
distribution of the microscopic configuration near each point x is a Bernoulli measure so 
that it is translation-invariant. 

(4) As we will see in the proof of Lemma \4.2[ a is increasing in p and b is increasing in m, 
respectively. In particular, the constant m measures the extent of the bias for the particles, 
that is, a larger m implies more particles on the right side. 

(5) The function /3(x) appears in other context: The asymmetric simple exclusion process 
on Z with jump rates p to the right and q to the left satisfying 0<p<l,p + q = l, has 
vpu a ,b) with b = logp/q as its invariant measure for every a G (0,1), where vpi.\ denotes 
the product measure on S such that £' I//3{ } [%] = (3(k),k G and the function j3(-;a,b) is 
defined by ((221) for all x G R {or x G Z), see [10], p. 382. 
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The macroscopic profile f3(x) has a connection to the Vershik curve which appears 
in a scaling limit for random Young diagrams distributed under the restricted uniform 
(Fermi) statistics. In fact, one can associate the height function tp e (u),u G [—£ — 1,£] of 
the Young diagram with the particle configuration n E by 

(2.8) ^(u)= 1*' 

fcSA^ :k>u 

Note that tjr is a right continuous non-increasing step function and satisfies 

(2.9) ^(-£-l) = K At ( V ), ^(£) = Q, 
with the area 

(2.10) / i,\u)du = {l+l)K kl {r,)+M kt (r 1 ). 

J-i-l 

Under the distribution v\^k,m, we consider the macroscopically scaled height function 
defined by 

(2.11) tf(x):=j^{lx), xG[-l,l]. 

Corollary 2.2. Under the conditions (|2.1|) and (|2.2[) . ifr converges as £ ^ oo to i)j in 
probability in the following sense: 



lim u A k,m sup V 0*0 - 4>(x) 
£—>oo 



xe[-i,i] 



>$ =0, 



/or every 5 > 0. T/ie Zimzt V * s defined by ip(x) = J^ 1 /3(y)dy,x G [— 1, 1] imt/t /3(x) 
determined in Theorem\2. 1[ 



From (gU, ([221), and (I2TTU|) . the limit ^ satisfies 



(2.12) ^(-l) = 2p, V(1) = 0, y ^(a;)dx = 2p + 4m 

and has a slope ip'(x) = —/3(x). Note that (|2.12p is consistent with (|2.5p . As we have seen 
in [3], Theorem 2.2, the Vershik curve iPr{x) defined for x > in the restricted uniform 
statistics having the area ipn(x)dx = 1 satisfies the ordinary differential equation 

(2.13) ^R + af/ R (l + ^ R ) = 0, 



where c = it /y/ 12. Here we consider in a rectangular box, and a simple computation 
with the help of (|2.6p shows that the limit in Corollary 12.21 satisfies the same ordinary 
differential equation (|2.13p with c = —b. Further discussions on the Vershik curves will be 
held in Section 5. 
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3 Local limit theorem for inhomogeneous Bernoulli sequence 
with unbounded weight 

For the proof of Theorem 12.11 we need to establish the local limit theorem jointly for the 
sum of inhomogeneous Bernoulli sequence and the sum with unbounded linearly growing 
weights k having some defects in them. Let a continuous function a = a(-) : [0, 1] —> (0, 1) 
and a sequence {a£ G (0, 1)}£ =1 satisfying the condition 

(3.1) lim max — a(k/n)\ = 0, 

n->oo l<fc<n 

be given. Let {Xk}^ =1 = {X'^}^ =1 be {0, l}-valued independent random variables with 
mean = for 1 < k < n and n G N. Note that a_ < a(x), < a+,x G 

[0, 1], 1 < k < n, n € N holds with some < a_ < a + < 1. We assume that a subset T n of 
{1,2,..., n} is given for each n G N and the size |T n | is uniformly bounded in n: \T n \ < C 
for all n. Under these settings, we consider the sums 

(3.2) S n = S n (X) := Y, X k, T n = T n (X) := kX k , 

for n G N, where = {1, 2, . . . , n} \ T n ; T n are defects in these sums. Then, it is easy to 
see that the (joint) central limit theorem holds for (S n ,T n ). Indeed, we define 

S n . — /tt~ (*~' ri En) > T n . — (T n Fji) ) 

with 

fcer= fcer- fcer= feer- 

where t>£ := a£(l — a£). Recalling that a(-) G C([0, 1]), the condition ([3.ip and |T n | < C, 
the asymptotic behaviors of E n , F n , U n and V n as n — > oo are given by 

(3.3) £ n = n(5 + o(l)), F n = n 2 (d + o(l)), 

C7 n = n(i; + o(l)), K = n 3 (7) + (1)), 

where a, a, v,v are positive constants defined by 

a = a(x)dx, a = xa(x)dx, 
Jo Jo 

.1 ,-1 

a(x)(l — a(x))dx, v= x 2 a(x)(l — a(x))dx, 



respectively. Then (S n ,T n ) weakly converges to (Yi, Y2) as n — >• 00, where Y = (Yi, Y2) is 
an R 2 -valued Gaussian random variable with mean and -E^Y] 2 ] = -EfYg 2 ] = 1, E\Y{Y2\ = A, 
where 

1 f 1 

X = —= / xa(x)(l — a(x))dx. 



'VV Jo 
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Indeed, the convergence of the corresponding characteristic functions is shown by Lemma 
3.21 below. Note that A G (0, 1) by Schwarz's inequality. The joint distribution density 
function of Y is given by 



(3.4) 



Ay) 



i 



2lTy/l - A 2 



exp 



y\ - 2X yi y 2 + y\ 
2(1 -A 2 ) 



, 2/ = (2/1,2/2) e 



We now state the corresponding local limit theorem. The set of all possible values of 
(S n ,T n ) is denoted by 

V n := {(K, L) G Z + x Z + ;P(S n = K,T n = L) > 0}, 

where Z+ = {0, 1, 2, . . .}. 

Proposition 3.1. We have that 



where 
(3.5) 



lim sup 

n ^°° (K,L)£p n 



2/1 



^/U^ n P{S n = K,T n = L) — q ( yi , 



(K — E n 



JJ2 



:(L-F n ). 



Remark 3.1. If a(-) G C 1 ([0, 1]) and ct^ are given by = a(k/n), then the convergence 
in Proposition \3. 1\ takes place with speed 0{\/y/n): 



sup 



^/U^y n P{S n = K,T n = L) 



) (2/1,2/2} 



< 



c 



n 



with some C > 0. However, in Section 4, we are forced to consider more general 
satisfying the condition (|3.ip . See Remark \4-l\ below. 



The local limit theorem for T n was studied by [8] in homogeneous Bernoulli case 
without defects and we extend it to the joint variables (S n ,T n ) in inhomogeneous case 
with defects. The rest of this section is devoted to the proof of Proposition 13. 11 We 
essentially follow the arguments in [11], but because of the unboundedness of the weights 
k for T n , a phenomenon different from the classical situation arises in the Fourier mode 
especially for the term I2 introduced below. 

Let f(s, t; S n , T n ), s, t G R, be the characteristic function of M 2 -valued random variable 

(Sri) Tn)' 



f(s,t;S n ,T n ) = E 



J(sS n +tT n ) 



where i 



-1. Then, for all (K, L) G V n , we have that 



(27r) 2 P(5 n = K, T n = L) = f f e-^ K+tL ^f( S , t; S n ,T n ) dsdt 

J —TV J —TV 



^UnTV 



V„1T 



-i(sy 1 +ty 2 ) 



f(s,t;S n ,T n )dsdt, 
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where y\ , y 2 are defined by (|3.5p . The second equality is due to a change of variables 
noting that 

(3.6) f(s,t;S n ,T n ) = exp{i(sE n + tF n )} f (y/TI^s , y/Vnt; S n ,f n ). 

Thus, if we define the error term by 



R n (K,L) = (27T) 2 { V r U^P(S n = K,T n = L) - q ( yi ,y 2 )}, 
it can be rewritten as 

/yfOnir r\/V^ir 
/ e-^+^fis^-Sn^dsdt 

e -i{s Vl +ty 2 ) & -{s 2 +2\st+t 2 )/2 dsdt 



/R 2 

since (27r) 2 (/o(yi> 2/2) is given by the second integral in the above expression. We divide 
R n (K,L) into the sum of three integrals: 

R n (K,L) = h + I 2 + h, 

where 

h = [ e-'^+to) \f( s ,t; S n ,f n ) ~ e -^ +2Xst+t ^ 2 } dsdt, 
h= I e- i( -^ +t ^f(s,t;S n ,f n )dsdt, 

JD 2 ,n 

h = I e-^+^e-^+^+^^dsdt, 
JD 3 , n 

respectively. Three domains are defined by 

Di, n = {{s,t) £ M. 2 ;\s\,\t\ < c^}, 

F>2,n = {(•?, t) G M 2 ; cy/n < \ s\ < \/Unir or cyfn < \t\ < \/Vr l n}, 
Dz, n = {( s >*) 6 l s l — c\fn or \t\ > cy/n}, 

respectively, with a small enough c G (0, \/^tt) chosen later. 

The estimate on ^3 = I^^ n {K,L) is easy. In fact, noting that s 2 + 2Xst + t 2 > 
(1 - \)(s 2 + t 2 ) and A < 1, we have a uniform bound on 23: There exist C\,c\ > such 
that 

(3.7) sup|/ 3 | < C ie ~ cin . 

K,L 

To give the estimate on I±, we prepare a lemma which is a two-dimensional version 
of Lemma 1 in Chapter V of p. 109. 

Lemma 3.2. For every 5 > 0, there exist no £ N, c 2 , C3 > suc/i i/iai, if n > uq, 



f(s,t;S n ,T n ) - e -(* 2 +2Asi+i 2 )/2 < + j_ s , + ,, )f 

ZioWs /or ewery (s,i) G M 2 satisfying |s|, |i| < c^y/n. 



|3 , c 2 , ,2n -c 2 (s 2 +t 2 ) 
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Proof. A simple computation leads to 



(3.8) 



f(s,t; S n ,T n ) = Yl i n /, r 

fcer: 



«7fe(l-o/=) 



+ (l-a fc )e-^ Qfc }, 



where 7& = 7^(s,t) := s / y/Un + tk / y/Vn and we simply write a k instead of ajj. If |i| < 
C3v^ with C3 > chosen later, from (|3.3p . 7^ can be estimated as 

(3.9) |7 fc | < cc 3 , 

with some c > 0. Therefore, since Taylor's formula implies 

1 



1 + z + -z z 



< c 4 i^r, 



with some C4 > for all z 6 C : |z| < CC3, we have that 



(3.10) a k e 



»7fc(l-«fc) 



+ (1 - a fc )e" 



= a fc (l + «7fc(l - «fc) - -tI(1 - ttfe) 2 ) + (1 - «fc)(l - ilk&k ~ ^lk a k) + R k 

1 Vk 2 , r> 
= 1 - y7 fc + #fc, 

with an error term R k G C having an estimate: 

(3.11) |22 fc | < C 4 (a fc (l - a fc ) 3 |7fc| 3 + (1 - a k )a 3 k \ 7k \ 3 ) < C 5 \ 7k \ 3 , 

where we write v k for v%. Let log(l + z) be the principal value defined for z £ C \ {z = 
iGR;i< — 1}. Then, 

|log(l + *)-z| <C 6 |z| 2 
holds for \z\ < 1/2 with some C§ > 0, and therefore, from (|3.8j) and (|3.10p . we have that 



(3.12) 



log f(s,t;S n ,f n )- (-y7fc + #fc) 



fcers 



< 



if I - v k ~/l/2 + R k \ < 1/2 for all 1 < k < n. Note that, from and (ETTTI) . the last 
condition holds if we choose C3 > sufficiently small. However, we see that 



(3.13) 



v kl 2 = (s 2 + 2Xst + t 2 )(l + o(l)), 

fcers 



by using (|3.3p and recalling |r„| < C. We thus obtain that 

1 



(3.14) 

with an error term R having a bound: 



log f(s,t;S n ,T n ) = -^{s 2 + 2\st + t 2 ) + R, 



n / 1 n n \ 

\r\ < m + 2c 6 1- J2 ^i7 fc i 4 + E i^i 2 + + * 2 )- 

fc=i V fe=i fe=i / 



The sum over is bounded by that over all 1 < k < n. However, from f)3.11|) and (j3.3|) . 
we have that 

n n „ 

(3-15) J2 W hk\ 3 < ~^(\s\ 3 + \t\ 3 ), 

fc=l fc=l ^ n 

with some C7 > 0. For the sum of f 2 .^! 4 , since v\ < (1/4) 2 , we have that 

n n ,4 i,|4,4\ 

(3.16) E^i 4 <^E(^ + ^ L ) 

k=l k=l v n n J 

<^(M 4 + ltl 4 )<^(M 3 + l*l 3 ), 

n vn 



with some C%, Cg > 0, if |s|, |i| < 03-^/71. Since ()3.9[) and ()3,lip show that is bounded, 
(|3,15p implies that 

n „ 

(3-17) 5> fe | 2 <%(| S | 3 + |f). 

fc=l v 

Therefore, (|3.15p - (|3.17p can be summarized into 

\R\ < ^(\s\ 3 + \t\ 3 ) + o(l)(s 2 + t 2 ), 



n 



if \s\, \t\ < csy/n. Coming back to (|3.14p . we have that 



/(s, t; S* n , T n ) e 



-(s 2 +2\st+t 2 )/2 



_ e -(s 2 +2Ast+t 2 )/2| e R _ 1 | 
< e -( 1 - A )( s2 + i2 )/ 2 | J R|| e l fl l. 



However, if |s|, \t\ < c^y/n, \R\ < (C11C3 + o(l))(s 2 + t 2 ) and therefore, by choosing C3 > 
sufficiently small and no sufficiently large, we have that 

e -(i-A)( s Ht 2 )/2 e |tf| < Cl2 e- c ^ 2+t2 \ 

with some C2, C12 > for every n > uq. We have thus completed the proof of the lemma 
by changing the choice of no to bound \R\ by means of the given 6 > 0, if necessary. □ 

Lemma 13.21 gives a uniform estimate on I\ = Ii tn (K,L) under the choice of c = C3: 
For every 5 > 0, 

(3.18) snp\h\<8 [ (\s\ 3 + \t\ 3 + s 2 + t 2 )e- C2is2+t2) dsdt, 

K,L 7R2 

holds if n > no; note that the last integral is converging. 

Finally, we give an estimate on I2 = l2 n (K,L). Using the relation (|3,6p . I2 can be 
rewritten and estimated as 

(3.19) \h\<\/UnVnf \f{s,t; S n ,T n )\ dsdt, 

JE 2 , n 
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where E% >n = {(s,i);c4 < |s| < ir or c^/n < \t\ < ir} and C4 = c inf„ e N y n/U n , 
C5 = c inf ng N > 0. Once we can show that there exist sufficiently small 9 > and 

k G (0, 1) such that 

(3.20) tt{A;;l < k < n,k $ T n ,s + kt $ [-9,9] mod 2tt} > ku - 1 - |r n |, 

for all (s,t) G Ez m then we have that 



\f{s, t; S n , T n ) 



IT { ak ' 



< r 



m- 1— C 



^ s+fe *) + (1 - a k )} 

where = and 

by recalling < a_ < a& < a + < 1. This together with (|3.19p and (|3.3p proves 



r := max max \a k e LU + (1 — < 1, 
M:0<|«|<7rfcer% 



(3.21) 



sup I h I < (2vr)V^nKr K ' 
K,L 

< C'n 2 r Kn - l - c . 



nn—l—C 



Three uniform estimates (|3.7|) . (|3.18|) and (|3.2ip conclude the proof of Proposition 13. 1[ 

The final task is to show (|3.2Up for all (s,t) G ^.n- To this end, we may assume 
t > by symmetry. We may also assume T n = and prove (|3.20p without — |r n | in 
the right hand side. We rewrite the region E2 n H {t > 0} into a union of three regions: 

E 2 , n n {t > 0} = #g u u 41 where 



2,n 
(2) 



2 , n -{(«,*); ^ <i<20,|s| <vr}, 



(3) 
2,n 



{(s,t);0 < t < — ,04 < |s| < vr}. 



7? 



Note that c(= C3) and therefore C5 was chosen sufficiently small so that we may assume 
< c 5 < 2vr. For (s,t) G since + G [-0,0] mod 2vr" implies "s+(k+l)t £ [-9,9] 

mod 27r", it is obvious that (|3.20p holds for T n = and k = 1/2. Now we take (s, t) G E\ n - 
The n points {kt}™ =1 are arranged on R in an equal distance and the interval [t, nt] 
containing all these points are covered by at most m := [(n— l)i/(27r)]+l disjoint intervals 
of length 2ir, where [ ] means the integer part. However, for an arbitrary interval Jet 
of length 27T, 

26 

${k; l<k<n,ktel,s + kt€ [-9, 9] mod 2vr} < — + 1. 

(2) 

Thus, since c^/n < t < 29 for (s, t) G E% „, we have that 

U{Jfe; 1 < k < n, s + kt G [-9, 9] mod 2vr} < m f y + 1^ < + lj + 1^ 
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by choosing 9 : < 9 < (1/vr + l/c 5 )~V 4 ) and this proves (pT20D for r n = and k = 1/2. 
Finally we take (s,t) G E^. Then, since < C5 < 2ir and < t < c§/n, n points {kt}^ =l 
are all located in the interval [0, 2tt). Now choose 9 : < 9 < C4/8. Then, recalling that 
C4 < \s\ < 7r, for example in the case that — 7r < s < — C4, 

{fc;l<Kn,s + tte [-0, 6>] mod 2vr} = {A;; 1 < k < n, s + kt G [-0, 0]} 

C {A;; 1 < k < n, kt > -9 - s} C {k; 1 < k < n, k > —n}. 

8c 5 

Thus we obtain (|3,20p for T n = by choosing k = (704/(805)) A 1 > 0. The case that 
C4 < s < ir can be discussed in a similar way. The proof of ()3.20p is completed for all 
(s,t) G E 2 , n - 

Remark 3.2. Our argument relies on the specific form k of the weights. It can be extended 
to linearly growing weights, but not for general weights such as a power of k. 



4 Proof of Theorem 12.1 

For a function /3 = /?(•) : [—1,1] — > (0,1), we denote by v^f-. the distribution on £a £ of 
{0, l}-valued independent sequences {r]k}keA e such that E[r]k] = /3(k/£), k G A^. The next 
lemma explains the reason that the functions of the form (|2.4p appear in the limit. 

Lemma 4.1. For a function /?(•) = /3(-;a, b) of the form (|2.4p . the conditional measure 
of v pf\ 071 ^A e ,K,M is a uniform probability measure for every a G (0, 1) and b G M: 

u ^[.)('\^A£,K,m) = u A e ,K,M(')- 

Proof. For a G (0, 1), let [i a be the probability measure on {0, 1} defined by /u a (l) = a. 
Then, if a = e b a/ (e b a + (1 — a)) for some a G (0, 1) and b G R, it holds that 

Ma(£)=^ 6 VW£) 5 

for £ = 0, 1 with a constant z a ^ = e b a + (1 — a). Therefore, for Vpf\ on £a^ with /?(•) of 
the form (|2.4p . we have that 

= n = n *i// few w%) 

= Z -l e feE fe6A ,^ fc /^A, (r?)) 

for 7] G Sa^ with a normalizing constant Z = HfeeA z abk/e = E A e [e b '^-' keA t kr,k ^ e ]. This 
implies that 

u P(.)(.'\^A-i>K,m) = ^("I^a^m), 
since YlkeAt^Vk = on Sa^^.m- However, u^ l {-\K\ t = K) is a uniform measure on 
{?7 G Ea^-Ka^) = so that Va l {'\^A t ,K,M) is also a uniform measure on T>A e ,K,M, that 
is, VA t ,K,M- □ 
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We next establish the one-to-one correspondence between (p, m) and (a, b) defined by 
(|2.5|) , In particular, one can uniquely find (a, 6) for every (p, m) in Theorem 12, li Consider 
the map $ for (a, b) G (0,1) x I to (p,m) = (F(a, 6), G(a, 6)) & D := {(p,m);p G 
(0, 1), m G (— u/2, u/2)} with u = p(l — p) defined by 



Lemma 4.2. The map ^ is a diffeomorphism from (0, 1) X R onto the domain D. 

Proof. Recalling that f3(x; a, b) is given by j5{x; a, 6) = e bx a/g(x; a, b) with g(x) = g(x; a, 6) = 
e bx a + (1 — a), we can easily compute the derivatives of F and G: 



This, with the help of Schwarz's inequality, implies that the Jacobian of the map <3? is 
positive everywhere, that is, J = — > 0- I n particular, the map $ is a local 

diffeomorphism. 

For every b £ R, set C b := {(F(a,b),G(a,b)) £ D;a £ (0,1)}. Then, C b is a Jordan 
arc in D connecting two points (0,0) and (1,0). Indeed, from dF/da > 0, Cj has no 
double points so that it is a Jordan arc. As o | 0, we have /3(x; a, b) — > so that 
(F(a,6),G(a,6)) -> (0,0), while (F(a,b),G(a,b)) -> (1,0) as a f 1, since 0(x;a,b) -»■ 1. 
Especially, Co = {(a, 0);a G (0,1)} is a line segment connecting (0,0) and (1,0). From 
the fact that J > 0, dF/da > and dG/db > 0, we see that the arc C bl is located above 
Cb 2 in the p-m plane if 6i > 62. This proves that the map $ is one to one. 

To show the onto property of we consider (a, b) satisfying F(a, b) = p for a fixed 
p G (0, 1). From (|2.7p . this condition can be rewritten as 




/3(x; a, b)dx 





e 2b P _ 1 



a 



e & _|_ e 2bp _ g 6(2p-l) _ ^ 



and therefore 



/3(x; a, b) 



e bx( e 2bp _ y 



,bxtp2bp 



1) + (e b - e^P- 1 ) ) 



which behaves as 




x + 2p > 1 
x + 2p < 1 



as 6 — > 00. Thus we have 




under the condition F(a, b) = p. In a similar way, we can show that 



lim G(a,b) = --p(l-p). 

— >— 00 z 



This proves the onto property of the map 



□ 
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We are now at the position to complete the proof of Theorem 12.11 Our method is 
standard in the sense that we apply the local limit theorem to compute the conditional 
distributions, see, e.g., [9], p. 353. The novelty lies in the following point. When applying 
Proposition 13.11 the term <Zo(2/i>2/2) needs to be uniformly positive to dominate the error 
term and this restricts our choice of the underlying measures in such a way that both y\ 
and ?/2 are sufficiently close to 0. For instance, in [9], it was only required to adjust the 
density parameter so that the macroscopic profile was constant over the space, but here 
the problem involves two parameters K and M or p and m. However, once we take the 
inhomogeneous Bernoulli measures with the functions f3(-;a,b) as macroscopic profiles, 
our choice allows two parameters a and b. This suits our purpose and we can realize the 
situation that both y\ and yi are close to at the same time under a suitable choice. The 
local limit theorem with defects is prepared to treat the numerator Bi introduced later. 

We denote the supports of the local functions fj in Theorem 12.11 by Fj C Z and set 
r(= rw) = U^jTfc.Xj, where T^Tj = Tj + kj and recall kj = kjj. Note that rcA< and 
{t^T j} P j = \ are disjoint if t is sufficiently large, since kj asymptotically behave as Xj£ and 
{xj}-j =1 C (—1,1) are distinct. Then, for every function /?£(•) of the form (|2.4p . Lemma 
14.11 shows that 



(4.1) E UWI J] r k] fj] - J] E mxj) [f 3i 



^(t/If = Z,K Ae ( V ) = K,M Ae ( V ) = M) 
At 



~ (11 /?(£?) 11 / '": . : Jj) \, , 

€£{0,1}^ j=1 j=1 v m\ K ^M = K,M Ae {rj) = M) 

where /3(-) is the function determined in Theorem 12.11 t]\y stands for the restriction of r/ 
to T and we denote by £j = £|r-. For given K = Ki and M = Mi, our special choice of 
j3i{-) will be the function /3^(-) = f3(-;ae,be) of the form (|2,4p with the solution (a^,^) of 
two equations <^E) for p = K/(2£ + 1) and m = M/{2£ + l) 2 . 

To apply Proposition 13. II for {r/k}k£A e , we need to shift it and consider X = {Xk}^ =1 
with n = 21 + 1 determined by = %-^-i(= {r^^rfjk) for k = 1, 2, . . . , n. For such X 
and defects T n C {1, 2, . . . , n}, we denote two sums S n and T n in ([32]) by S^ n (X) and 
T^ n (X), respectively, to indicate the defects clearly. Then, we have that 

(4.2) S 9 n (X) = K Ae ( V ) and T n (X) = M Al {r,) + {£ + l)K Ae ( V ). 

The numerator of the fractional expression in the right hand side of (|4.ip is equal to 

P a n (X\ Te+lV = r-\e, S*(X) = K, T®(X) = M + (£ + l)K) 



r (OPan (s T n e+ir (X) = K — st +irr (r- + \0, 



Tl t+lT {X) = M + {£ + 1)K — T^ +ir)C (r- + \e) 

where a n = {a^}^ =1 := {fii{{k — £ — l)/£)} k=1 - We have denoted by P a n the distribution 
of X such that E[Xk] = for 1 < k < n, by ^(.) the measure ^5.) restricted on T, and 
{ti + \T) c = {1, 2, . . . , n} \ T£ + iT. However, under the scaling conditions (|2.ip and ()2.2p . by 
Lemma |4.2| /3^(-) converges to /3(-) uniformly: 

(4.3) lim max W e (x) - /3(x)| = 0. 

e^oo xe{-i,i] 
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Hence, we have that 



E - n^)[/;K(#) = «> 

?e{o,i} r i =1 i =1 



as ^ — > oo. Therefore, (|4.1|) can be rewritten as 

p v 



(4.4) 



where 



E - Vi^MYm^) x ^ - 1 [+0(1), 



$e{o,i} r j= 



i=i 



^ := P a n [SI(X) = K, Tl{X) = M+(£+l)K 
Be := P a n (s T n ^ T (X) = K — st +lT) \r^), 



T^ V (X) = M + (£ + l)K - T^ +ir)C (r- + \0) . 



We now apply Proposition 13.11 to compute the asymptotic behaviors of At and B^. 
Note that, from (|4.3p . a n = {a&}fc =1 satisfies the condition (|3.ip with a(x) := (3(2x — 
1), x G [0, 1]. By the choice of (a^, 6^), we can show that yi = V2 = 0(1/ \/n) as £ — >• oo (or 
equivalently n — )• oo) for y\ and y2 defined by (|3,5p for both ^ and Bg. Indeed, for Ai, 



1 



yi 



1 



- E n 



Un I fe=l 



1 



/v - ([ I p t (x)dx + (J 



if 



2 * +0 a 



2£ + 1 



o(i/v^). 



Here, for the third equality, we have used the uniform bound: sup^ max rEg [_ 11 ] |/3^(x)| < oo, 
recall the choice of f3i(-) for the fourth and (|3.3p for the last. Similarly, 



V2 



fVr. 
1 



■■(L-F n ) 



1 f n 

Vn i fc=i 



k-£- 1 



M + (£ + 1)# - ^ 2 |y (x + l)ft(x)da; + O Q 



= 0(l/v^). 

For compared with Ai, there are additional terms 



1 'st^irr 1 



- E ' 

feer f+1 r 



15 



in y\ and 



ftio- E 



k-l-1 

£ 



in ?/2) but both these terms behave as 0{l/y/n). Thus, from Proposition 13.11 two terms 
Ai and Bg both behave as 



Ap, Bp 



1 



+o(i: 



as ^ — > oo (or equivalently n — >■ oo). Since 1 — A 2 > 0, this shows that Ai/Bg — >■ 1 and 
completes the poof of Theorem 12.11 from (|4.4p . 

Remark 4.1. 7/u>e ta&e /?(•) itself instead of (3t{), the terms y\ and 1/2 in the exponential 
of qo behave as o(y/n) so that qo in general converges to and the local limit theorem turns 
out to be useless. 



5 Proof of Corollary 12.21 and relations to Vershik curves 

For the proof of Corollary 12.21 it suffices to show the weaker convergence 



(5.1) 



lim vk,,k,m 



> 5 



0, 



for every 5 > and ip € C([— 1, 1]), where (ip, ip) = J_ 1 ip(x)(p(x)dx. In fact, (|5.ip implies 
the stronger convergence result stated in Corollary 12.21 due to the monotonicity of ijr, see 
Remark 2.5 in [3]. For proving (|5.ip . it is enough to show that (t/S , <p) converges to (ip,<p) 
in L 2 -sense. However, recalling ()2.1ip and (|2.8p . a simple computation leads to 



<PM = \ E 

k=-t- 

where (p(x) = J* (p(y)dy. Therefore, once we can show that 

1 y> fk 

1 2^ 



(5.2) 
and 
(5.3) 



lim E v 



Vk<P 



k=-e+i 



lim E u . M 



. k=-e+i 



' k 



</3,93> 2 , 



the L 2 -convergence follows by noting that (/3, (p) = (tp, 99). We only give the proof of (|5.3p . 
since (15.21) is similar and easier. To this end, from the estimate 



<2 sup \<p(xi) - tp(x 2 )\, 

\xi-X2\<j 



(ij; ,tp)- I r] [ex]+1 ip(x)dx 
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it is enough to prove 



(5-4) lim E UAi KM 



l— >oo 



as 

fl rl 



-\2 



^V[£x]<f( x ) dx 

where [£x] stands for the integer part of Ix. However, the expectation in (|5.4p is expanded 



'-1 J-l 

and Theorem 12.11 applied for p = 2 and f\ {rj) = f 2 (r/) = t]q implies that 



Thus, Lebesgue's convergence theorem shows (|5.4p and this completes the proof of Corol- 
lary H2J 

We finally compare our result in Corollary 12.21 with those in pQ . In [1] , the grand 
canonical ensembles in a rectangular box for the uniform (Bose) statistics are dealt based 
on a combinatorial method, while we have discussed the canonical ensembles in a rectan- 
gular box for the restricted uniform (Fermi) statistics due to a probabilistic approach. 

Theorem 1 of [T] shows, by rotating the plane coordinates by 45 degree, that the limit 
curve t G [0, 1] i— > L(t) in the box with the ratio of height/width given by p/(l — p), p £ 
(0, 1) is determined by 

msLfAt) = i - log m r 

where h(t) = e~ ct — e ct + e ~ c ( 2_2 ^~*) — e~ c ^~ 2 P\t G [0, 1], and c G R is a parameter which 
controls the area. By rotating back to the original coordinates, this curve is transformed 
to the curve v = 4>(u) in u-v plane given implicitly by 

(5.5) u = ^=(t + L(t))e[0,(l-p)/V2], v = -L(-t + L(t))e[-V2p,0}. 

However, as seen in Proposition 4.4 of [I], the curve v = <j){u) appearing in the uniform 
statistics can be related to the curve y = ip(x) in x-y plane appearing in the restricted 
uniform statistics by 

fit \- T/( \ -^(s)) 
p(x) — —w (X) = j , 

where G^(u) = u—<f>(u) and G^ 1 is its inverse function. But, (|5.5|) shows that G^,(u) = \^2t 
and thus, setting x = G^(u) G [0, y/2], we have that 

or equivalently 
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A simple computation leads to L'(t) = h! (t) / (ch(t)) and L"(t) = c(l — L'(t) 2 ), since 
h"(t) = c 2 h(t). This proves that 

Pi?) = --±=L"(x/V2) = ^= (L'(x/V2? - l) = -y/2cp(x)(l - Rx)l 

so that we arrive at the ordinary differential equation for ip: 

(5.6) 4>"(x) + v/2c^'(x)(l + i?(x)) = 0, xe[Q,V2]. 

Moreover, noting that L(0) = and L(l) = 1 — 2p, the height difference of ip at two 
boundary points is given by 

0(O)-V(\ / 2)= / P(x)dx = V2p, 



or, by normalizing ip(y/2) = 0, we have that ^(0) = \[2p. 

To compare ip with ip in Corollary 12.21 note that ip is defined for x £ [—1,1]. The 
shift in x does not change the form of the equation. To consider on the interval of same 
length, we introduce the scaling for ip defined on [0, y/2] by 

$r(x) : = l^( 7 s), xG[0,V2/ 7 ], 
7 

for 7 > 0. Then, if ip satisfies ip" + Kip'(l + 0') = 0, ip 1 satisfies the equation (ip' 7 )" + 
^(ip 1 ) 1 (1 + (ip 1 )') = 0. Applying this for the equation (|5.6|) with 7 = l/v2 and k = y2c, 
we can derive the equation for ip 1 /^ 2 : 

(^VV5)» + c(^/^)'(i + (^ 1/V5 )') = 0, i£ [0, 2], 

with ^/^(O) = 2p and , 1 / v/ ^(2) = 0. This coincides with the ordinary differential 
equation (|2.13|) for the Vershik curve with c = c and p = p. Thus, we can identify the 
limit curves in a rectangular box for grand canonical ensembles in the uniform statistics 
and for canonical ensembles in the restricted uniform statistics, namely ip 1 /^ 2 = ip (except 
the shift in x by 1), if the relations c = —b and p = p hold. 
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